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Spectral Index as a Function of Mass Accretion Rate in Black 
Hole Sources. Monte-Carlo Simulations and an Analytical 

Description 

Philippe Laurent^ and Lev Titarchuk^''^'^ 
ABSTRACT 

We present herein a theoretical study of correlations between spectral indexes 
of X-ray emergent spectra and mass accretion rate (rh) in black hole (BH) sources, 
which provide a definitive signature for BHs. It has been firmly established, using 
the Rossi X-ray Timing Explorer (RXTE) in numerous BH observations during 
hard-soft state spectral evolution, that the photon index of X-ray spectra in- 
creases when m increases and, moreover, the index saturates at high values of 
rh. In this Paper, we present theoretical arguments that the observationally 
established index saturation effect vs mass accretion rate is a signature of the 
bulk (converging) flow onto the black hole. Also, we demonstrate that the index 
saturation value depends on the plasma temperature of converging flow. We self- 
consistently calculate the Compton cloud (CC) plasma temperature as a function 
of mass accretion rate using the energy balance between energy dissipation and 
Compton cooling. We explain the observable phenomenon, index-m correlations 
using a Monte-Carlo simulation of radiative processes in the innermost part (CC) 
of a BH source and we account for the Comptonization processes in the presence 
of thermal and bulk motions, as basic types of plasma motion. We show that, 
when m increases, BH sources evolve to high and very soft states (HSS and VSS, 
respectively), in which the strong blackbody-like and steep power-law compo- 
nents are formed in the resulting X-ray spectrum. The simultaneous detections 
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of these two components strongly depends on sensitivity of high energy instru- 
ments, given that the relative contribution of the hard power-law tail in the 
resulting VSS spectrum can be very low, which is why, to date RXTE observa- 
tions of the VSS X-ray spectrum has been characterized by the presence of the 
strong BB-like component only. We also predict specific patterns for high-energy 
efold (cutoff) energy (Efoid) evolution with m for thermal and dynamical (bulk) 
Comptonization cases. For the former case, Efdd monotonically decreases with 
m, in the latter case, the E'/oi^— decrease is followed by its increase at high values 
of rh. The observational evolution of Efoid vs m can be one more test for the 
presence of a converging flow effect in the formation of the resulting spectra in 
the close vicinity of BHs. 

Subject headings: black hole physics — accretion disks — radiation mechanisms: 
nonthermal — X-rays: general 



Introduction 



A study of the characteristic changes in spectral and variability properties of X-ray 
binaries has proven to be a valuable source of information on the Physics governing the 
accretion processes and on the fundamental parameters of black hole (BH) sources. BH 
observational appearance is conventionally described in terms of BH state class ification (see 
Remillard fc McClintockl bood : iBellonil Eoosi iKlein-Wolt fc van der KlisI boosl . for different 
definitions of BH states). We adopt the following BH state classification for five major BH 
states: the quiescent, low-hard (LHS), intermediate (IS), high-soft (HSS) and very soft states 
(VSS). When a BH transient goes into outburst, it leaves the quiescent state and enters 
the LHS, a low luminosity state with the energy spectrum dominated by a Comptonization 
component combined (convolved) with a weak thermal component. The photon spectrum in 
the LHS is presumably a result of Comptonization (upscattering) of soft photons, originat- 
ing from a relatively weak a ccretion disk, off electrons of the hot ambient plasma [see e.g. 
Sunyaev &: Titarchuk I (Il980[ ). hereafter ST80]. The HSS photon spectrum is characterized 
by a prominent thermal component which is probably a signature of a strong emission com- 
ing from a geometrically thin accretion disk. A weak power-law component is also present 
at the level of not more than 20% of the total source fiux. 



Timing an d spectral propert i es of an accreting BH are tightly correlated for a number of 



BH sources [see 



Vignarca et al. 



(12 0031) ari d a comprehensive study of t hese correlations was 



done by lShaposhnikov fc Titarchuk I (l2009[ ). lTitarchuk fc Seifina I (|2009[ ). hereafter ST09 and 
TS09, respectively]. Correlations between the spectral hardness (photon index) and the char- 
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acteristic frequencies of the quasi-periodic oscillations (QPOs) observe d in the lightcurves of 



BH s ources has been pr oposed to use as a tool to determin e a BH mass (ITitarchuk. &: Fiorito 



2004 hereafter TF04). iRemillard fc McClintock I fl2006f ) demonstrated and then ST09 con- 
firmed that there is a correlation between values QPO frequencies and the disk flux, namely 
mass accre tion rate M in the disk, for a number of BH candidates (BHCs). Moreover, ST09, 
TS09 and ISeifina &: Titarchuk I (|2010| ). hereafter STIO, also found strong correlations be- 
tween spectral index and m, where m = M / MecLcI is dimensionless mass accretion rate in 
units of critical mass accretion rate MEdd = ^Eddl^? ■, in 19 spectral transition episodes from 
10 BH sources observed with the Rossi X-ray Timing Explorer (RXTE). ST09 combine these 
index-QPO and index-m correlations to measure BH masses in 7 BH sources. For some of 
these sources, (i.e. H 1743-322, Cyg X-1, XTE 1650-500, GX 339-4 and XTE J1859-226) this 
new ST09 method provides a much better precision than conventional dynamical methods 
(see ST09). In addition TS09 and STIO using index-m correlations evaluate BH masses in 
GRS 1915+105 and SS 433 respectively. 

The index-QPO and index-m correlation has a specific shape characterized by a rise 
part followed by the saturation plateau. The ST09 method is based on scaling (sliding) the 
correlation pattern related to a reference source vs that of a given source for which a BH 
mass is determined. This determination can be only done for this pair of the sources if the 
saturation values of index and inclination of the rising part are the same for both. Thus, the 
correlation patterns need to be self-similar to i mplement this scaling and, fi nally, to obtain a 
BH mass for a given source. Using results from lLaurent fc TitarchukI (jl999[ ). hereafter LT99, 
TF04 who were the first to suggest a model for the index-m correlation observed in BHCs, 
we further explore the model of index-m correlation and present a thorough modeling of this 
correlation using Monte Carlo simulations for a wide range of basic model parameters. We 
also investigate the possibility that the shape of the correlation pattern provides a direct 
signature of the bulk moti on (converging) flow onto a black h o le, wh i ch would be the signa- 



ture o f the black hole [see ITitarchuk. Mastichiadis fc Kylafis I (jl997| ) , ITitarchuk &: Zannias 



( 119981 ) and LT99 for more details on this subject]. We present theoretical arguments based 
on the Monte-Carlo simulation and analytic consideration that the index-mass accretion rate 
saturation effect observed in a number of BHCs is a signature of a BH converging flow. This 
saturation effect should occur when the mass accretion rate M exceeds the Eddington li mit 
AifEdd, which can only exist in BH sources, see ITitarchuk. Mastichiadis fc Kylafis I (119971 ). 



The origin of the spectral state transition and its fin al end as a steep power-law in th e 



high soft states is still debated in the literatu r e [see e .g. IChakrabarti fc Titarchuk I ( 



hereafter CT95, LT99, 



Laurent fc TitarchukI (120011 ). iNiedzwiecki fc Zdziarski 



995|) 



mm and 



Remillard fc McClintock I ( l2006h . hereafter RM06]. However, there is agreement within the 
community that the spectral transition is driven by mass accretion in the system (see e.g. 
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CT95, RM06) such that, when rh increases a BH source goes from the LHS to soft states (HSS 
or VSS) through the intermediate state (IS), but a detailed scenario for this transition is still 
missing. A natural question to ask is: why does the source go to either the HSS or the VSS 
and what physical processes are behind of each of these transitions? The soft component 
of the observed X-ray spectra is usually fitted by a blackbody (BB) shape modified by 
Comptonization for which the Comptonization parameter Ygf decreases, or the energy index 
of the Comptonization Green's function agj increases towards a softer state [see details in 
TS09]. Furthermore, when the mass accretion rate increases up to some critical value photon 
index Tg/ = a^z + l starts to saturate to a value of about 4.2. What mechanism is responsible 
for the behavior of Tgf vs rh and for the saturation level of the index? In contrast, ST09 
and TS09 show that the index of the hard component of the X-ray spectrum Th saturates 
to different levels that varies between 2 to 3 for different sources and for different outbursts 
in a given source. For example, saturates to 3 in GRS 1915-1-105 (TS09) and to 2.3 in SS 
433 (STIO) and it varies from 2 to 2.6 for different outbursts in GX 339-4 (see ST09). 

Another phenomenon revealed in X-ray observations of compact objects [neutron stars 
(NSs) and BHs] is related to the evolution of a high energy efold (cutoff) energy Efoid of 
the spectra during the spectral evolution. For NSs, Ef„iA s t eadily decreases from t h e har d 
state to the soft state [Titarchuk &: Shaposhnikov I ( 120051 ) . iFarinelli fc Titarchuk I ( 120101 )] 
but, in BHs, the observational pattern of Efoid vs m is such that Efoid begins to de- 
crease with rh, which is then follow ed by a shift of Ejom to higher values when rh increases 
( iTitarchuk fc Shaposhnikov 1120101 ). 

We present our picture of accretion in the innermost part of source, the Comptonization 
region [Compton cloud (CC)] and a model to calculate the CC plasma temperature in ^ 
We provide the details of the Monte-Carlo (MC) simulations in ^ We present a combined 
effect of thermal and bulk motion Comptonization and BH spectral evolution in ^ We 
show and discuss our MC simulation results and their theoretical explanation related to 
the index-mass accretion rate correlation in ^ Specifically, we discuss the signature of 
bulk motion Comptonization (BMC) and its relation to the index evolution during state 
transition. Also, we show that the index saturation effect is a direct consequence of the 
existence of the innermost bulk motion region and, therefore, can be considered to be an 
observational signature of a BH. Efoid-fh correlations results and their interpretation are 
shown in ^ Conclusions follow in ^ 
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2. The model 

We illustrate our accretion scenario in Figure [H which we suggest taking place in the 
innermost part of a black hole source. We assume that an accretion flow onto a black hole 
(BH) con sists of three parts: a g eome trically thin accretion disk [standard Shakura-Sunyaev 



disk, see IShakura fc Sunyaev I (Il973l )]. a transition layer (TL), which is an intermediate 
link between the accretion disk, and a converging (bulk) flow region (see TF04), that is 
assumed to exist, at least, below 3 Schwarzschild radii, 3rs = QGMhh/c^- There is substantial 
obs ervational evidence for this three componen t model for accretion flows onto black holes 



[see iTitarchuk. Shaposhnikov fc Areflev I (l2007f ). hereafter TSA07, ST09, TS09 and STIO]. 
In addition, power density spectra (PDSs) during the transition from the low-hard to high- 
soft states indicate low and high frequency white-red noise (WRN) components, with discrete 
Lorentzian QPO features superposed. The low frequency WRN component dominates in 
the high-soft state PDSs while the high frequency PDS component dominates in the low- 
hard state, with the interpretation that these two components represent the diffusion of 
perturbations in an outer accretion disk and an inner Compton corona (transition layer), 
respectively. These two regions are probably separated by an accretion shock, which moves 
to smaller radii in the high-soft state and to larger radii in the low-hard state (TSA07). This 
accretion shock represents a location in the disk where the outer Keplerian flow begins to 
transition to the inner sub-Keplerian motion of the compact object and its inner corona. 

We assume that the plasma temperature of the converging flow is near that of the TL 
(or corona). The TL temperature can be presented in the analytical form using Eq. (12) in 
TFO^E: ^ 

10^ " " ro(l + 3ro/8)(l + g)-0.62 

where tq is TL optical depth and q = Qd/QrL is a ratio of the energy releases in the disk 
and transition layer, respectively. 

We relate tq to optical depth of the converging flow Tff region by using the flow conti- 
nuity equation 

Vff 



where the free-fall velocity 



^ t/ \r=r.utrff, (2) 



Vfj = c{rs/rY'\ (3) 



c is the speed of light. Tout is an outer radius of the converging flow and < Vma > is a TL 
magneto-acoustic velocity averaged over the layer < Vma > and approximately equals to 



1 There is a typo in formula (13) of TF04 which is the solution of Eq. (12) in TF04. 
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radial velocity of the TL flow. 

The optical depth of the converging flow, as evaluated by an observer at rest, is given 

by 

Tff = m(7r/2 - arcsin x^/'^), (4) 

where m = M/M-^d is a dimensionless mass accretion rate in units of the critical mass 
accretion rate M^d = Leci/c^ and Xout = i^out/^s is the dimensionless outer radius of the 
converging flow in units of Schwarzsc hild radius rs- Note , if the dimensionless BH spin a is 



less than 0.8 then 2.7 < Xout < 3 [see iLandau fc Lifshitz I (11975 I )] 



Titarchuk. Bradshaw &: Wood I (l200ll ). hereafter TBWOl, estimated the i?-fleld around 
neutron stars (NS) based on observed kilohertz and viscous quasi-periodic oscillation (QPO) 
frequencies. They found that the best-flt values of Vma in the transition layer related to the 
i?- field of NS and plasma density is in the range of 10* — 10^ cm s~^, where 
Vma~^0^{B /lO^ G)/(47rp/10~^g cm~^)^/^ cm s~^ and p is TL plasma density which is of 
order of p = neTUp = 1.6 x 10~^(ne/10^^cm~^) g cm~^. Unfortunatel y, there is no estimate o f 



Vma for the innermost part of a BH that is based on observations. iTitarchuk et al. I (Il998l ) 
introduced the Reynolds number Re = VrR/z> (7 in their notation) where Vr, D are an 
average radial velocity, an average viscosity over a given configuration, respectively and R is 
a configuration scale. They demonstrated that the size of the transition layer (CC) between 
the fast rotating accretion disk and the relatively slow rotating central object (either BH or 



NS) s trongly depends on the Reynolds number Re. Moreover, ITitarchuk fc Shaposhnikov 



(120081 ) infer, by analyzing X-ray data for Cyg X-1, that the Reynolds number Re of accretion 
flow in the transition layer increases from 10 to 100 when the source evolves from the low/hard 
state to the high/soft state. This is possible when an increase of Vr (~ Vma) leads to the 
Re rise. 

Because m is proportional to Vr (and ~ Vma) indicates that Vma is also proportional 
to Tff (see Eq. IHfor linear relation between rh and Tff): 

< Vma >= Co + CiVfj (5) 
where Cq ~ 0.01c and Ci/Cq -0-10. 

The combination of Eqs. ([2]), ([3]) and ([5]) gives tq as a function of Tff and Xout- 



[1 + (C,/C„)r„] 



X 



1/2- 

out 



Equations (jl]) and (|6]) yield, for nonzero Ci, TL optical depth tq saturates to a constant 
value, TQ^sat = c/{Cix]l^f-) when Tff (or fn) increases. This saturation of tq leads to the 
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saturation of temperature Tg with m for a given q (see Eq. [T]). In this case, the temperature 
saturation value (see Eq. [1]) is 

Te,.ai = i"(ro,sat,g)10'K. (7) 



The thermal Componization index a^/ als o saturates when rn and Tp sa turate given that a 



is a specific function of tq and Tg [see e.g. ISunyaev &: Titarchuk I (|1980[ ). hereafter ST80]. In 



■■f 



Figure [21 we present tq and Te (in units of keV ) as functions of m calculated using Eqs. ([T]) 
and ([6]). These calculations are made for q = 0.01, but for all g ^ 1 the results are similar 
to those presented in Figure [2l It is clearly evident in Figure [2] that saturation of tq and Te 
occurs at high values of m. 

Using Eq. and for tq > 1, the product TeTg 4 x 10^/(1 + q) K. The thermal 
Comptonization parameter Ysf oc {kTe/meC^)TQ (ST80), which implies that a^/ ~ 
goes to a constant too. In §U we demonstrate that this saturation effect of the thermal 
Comptonization index is reproduced in our Monte Carlo simulations. It is also worthwhile 
to emphasize that the index saturation of the soft component of X-ray spectra is detected 
in RXTE observations of GRS 1915+105 (see TS09). 



3. The Monte-Carlo simulations 

The parameters of the simulations are g, Ci/Cq. We fix value of Cq ~ 0.01c, Xout = 3 and 
q = 0.01 and we change rh from 0.01 to 10. As it is seen from Eq. ([T]) the plasma temperature 
Tg of the corona region is not affected if the disk cooling factor q <^ 1. Therefore, for most 
of the simulated cases, we fix the q value at 0.01. Once q and Ci/Cq are set, we compute 
Tff, Tq, and Te using Eqs. @, and ([1]), respectively. 

The physical model for the Monte-Carlo (MC) simulations is the Comptonization of the 
soft disk photons in the converging flow of optical depth Tfj extended from to Tout = Srg, 
surrounded by a spherical shell of optical depth tq, corresponding to the transition layer. 
The MC simulated spectra are not sensitive to the geometrical size of the TL layer, the 
radius of which being set to lOrs, but is very sensitive to tq, Tjf and Tg. Thus, we have 
incorporated plasma free-fall accretion onto the central black hole in the converging flow 
region as is described in LT99. Additionally, we include the thermal motion of the electrons 
in the converging flow and transition layer, simulated for electron temperature Tg given by 
Eq. (dD. 

The seed X-ray photons were generated uniformly and isotropically at the surface of the 
border of the accretion disk, from lOrg to 12rs. The soft seed photon spectrum is a blackbody 
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spectrum with a temperature of 0.9 ke V, which is close to the color BB temperature found 
in the spectra of BHC sources [see e.g. iBorozdin et al. I (119991 )]. 



In our MC simulations, we followed the trajectory of each photon in the following way: 
first, a uniform deviate U is generated such that Tph = — In f/, Tph being the optical thickness 
a photon will travel in curved geometry, we use the Schwarzschild metric, before scattering 
off an electron. We then integrated the optical thickness At = naAl along the photon path 
up to Tph, taking into account the variation of the cloud number density n and of the cross 
section a with the radius. Here the cross section a is the Klein-Nishina cross section averaged 
over the local relativistic thermal electron distribution. The gravitational red (or blue) shift 
endured by the photon was also computed at each step of this integration. 

We simulate photon Compton scattering off an electron that includes the exact motion 
of the electron. We also take into account the case where, at the end of the integration, 
the comput ed optical thickness has r eached Tph and the photon has not left the considered 
volume [see iPozdnyakov et al I (119831 )]. In our procedure, we first compute the scattering 
electron momentum, and derive the scattered photon and electron characteristics from the 
Compton scattering kinetics. Then we check to determine if this event is consistent with the 
Compton scattering probabilities; if it is, then the event is kept; if it is not, another scattering 
photon is generated, and the process goes on until the event is accepted. This process was 
su ccessfully checked by compar ing its results (see LT99) with the analytical results presented 



m 



Titarchuk &: Zannias I (119981 ). Once the new energy and direction of the photon has been 



determined by the Compton kinetics, we track it in the same way as above until it makes 
another scattering, it escapes from the considered volume, or until it is "absorbed" by the 
black hole at its horizon. 



4. A combined effect of thermal and bulk motion Comptonization. BH 

spectral evolution 

For a BH of the dimensio nless spin a < 0.8 th e radiu s of t he marginally stable orbit in 
the disk Vmso is about 3rs [see iRuffini fc Wheeler I ( 1l97ll ) and iLandau fc Lifshitz I ( 1l975 I )]. 
Thus, below r^so and if the plasma temperature of the flow is essentially non- relativistic, 
the accretion flow is almost in a free-fall regime, i.e. advection dominated, where the main 
effect of the up-scattering of the disk soft photons is the dynamical (or bulk) Comptonization 
Jhaurent fc TitarchuklE999l ). 



In Figure [3l we illustrate the effect of bulk motion Comptonization on the spectral 
signature. The parameters for the simulated spectrum were rh = 5, g = 0.01, and Ci/Cq = 7, 
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resulting in and inferred electron temperature of the flow kTe = 4.3 keV (see Eqs. [H |U |6]). 
The dotted line corresponds to the simulated spectrum where the bulk effect of the converging 
flow is not taken into account and the solid line corresponds to the spectrum where the bulk 
velocity of the converging flow is included in the simulations. The effect of bulk (dynamical 
Comptonization) flow is clearly seen above 20-30 keV. The steep power-law component that 
emerges in the spectrum which is now extended up to 200 keV. 

All our simulated spectra can be represented by an additive XSPEC model consisting of 
the addition of two, so called "bulk motion" Comptonization (BMC) components: a BMC 
with a high energy cut-off [BMCl component) and a relatively soft BMC2 component with 
a high energy cut-off. Thus, our model for fitting is [bmcl * highecutl + bmc2 * highecut2). 
The BMC model describes the outgoing spectrum as a convolution of the input "seed" 
blackbody spectrum, with a normalization of A^tmc and color temperature is kT, with the 
Comptonization Green's function which is a broken power law of spectral index a (photon 
index T = a + 1). The resulting spectrum is characterized by the parameter log(A) related 
to the Comptonized fraction / as / = A/{1 + A) and spectral index a. 

In Figure m we show an example of the simulated spectra produced by two BMC compo- 
nents with high energy cutoff, specifically bmcl*highcutl +bmc2* highcut2. The spectrum is 
shown in units of keV (keV cm~^s~^keV~^) where absolute values are multiplied by an arbi- 
trary constant. The hard component is given in dash green and the soft component is given 
in dash red. The resulting best-fit spectrum is presented by the solid green curve. In our 
simulations the mass accretion rate in Eddington units, is m = 5, Ci/Cq = 5, q = 0.01 and 
the flow electron temperature is kTe = 2.4 keV. The characteristics (parameters) of bmcl and 
highcutl are spectral index at = 1-8 (Th = 2.8), log A = 0.26, cutoff energy Ecut = 100 keV, 
efold energy Efou = 500 keV and that for bmc2 and highcut2 are a^/ = 2.23 (F^j = 3.33), 
log A = 0.99, E,^t = 4.5 keV, Ef^id = 4.2 keV. 

The general picture of LHS-IS-HSS transition is illustrated in Figure [5] where we bring 
together simulated spectra related to the observational spectra of LHS, IS, HSS and VSS to 
demonstrate the BH spectral evolution from the low-hard to soft states. We reproduce the 
observational spectral evolution by varying the mass accretion rate rh from 0.05 to 5. Top 
and bottom panels correspond to Ci/Cq = 5 and Ci/Cq = 10, respectively. The different 
shapes of the simulated spectra are related to the different spectral states in the following 
manner: low m = 0.05 (black histogram), which the TL optical depth tq = 2.5, 2 and 
kTe = 40, 55 keV for Ci/Cq = 5 and 10, respectively (see Fig. [2]), is similar to the observed 
LHS spectrum (see ST09, TS09); increasing the accretion rate to m = 0.1 (blue histogram), 
the spectrum becomes softer (i.e. power law is steeper and exponential cutoff energy is shifted 
to lower energies than that for m = 0.05) and the Comptonization component dominates 
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the blackbody (BB) component, which is barely seen in 3-150 keV energy range is similar to 
the spectrum seen in the beginning of the IS state; increasing rh further produces a spectra 
(green histogram) characterized by a strong soft BB component and an extended steep power 
law, which is comparable to the observed HSS spectra. 

For moderate m = 0.5 (red histogram) and Ci/Cq = 5 the spectral shape is a blackbody- 
like but slightly modified by the thermal Comptonization because kTf, ~ 4 keV is very 
low for efficient thermal up-scattering and Tjf ~ 1.6, m = 0.8 is also low for the bulk 
motion Comptonization. A similar spectra is detected during the very soft state (VSS). 
However when mass accretion further increases to a value of m = 5 the extended high energy 
component is seen up to energies 300 keV as a result of the bulk motion Comptonization in 
the converging inflow, r < r^. 

For Ci/Cq = 10 and m > 1 the spectral shapes are strongly affected by thermal Comp- 
tonization because kTe ~ 15 keV and tq ~ 8 — 10 are relatively high [see Figs. [2], [5] (lower 
panel)]. One can also observe the effect of bulk motion Comptonization seen as an extension 
of the emergent spectrum to energies of 300 keV when mass accretion rate increases, see 
green curve for m = 5 and compare the similar effect for Ci/Cq = 5 (Fig. O lower and upper 
panels respectively). 

5. Indices of soft and hard components of the X-ray simulated spectra as a 

function of mass accretion rate 

We have made series of simulations with different values of the Ci / Cq ratio in order to 
study the evolution of Ci/Cq - index saturation. This evolution is shown on Figures EHU 
where the saturation limits for Tsat.s and Tsat,h are shown. Moreover, we show the observed 
pattern of index vs fn (see ST09 and TS09) can be reproduced in our simulations by the 
variation of only one parameter of the velocity profile, i.e. Ci/Cq. In Figures [MHl we show 
Vgf (red points) and (blue points) as a function of m for a given Ci/Cq. Saturation 
of Tgf occurs in all simulated cases. The highest level saturation of F^j is achieved at 
Ci/Cq = 5, when the matter accumulation is strong in the transition layer (see Fig. [9]). The 
saturation value Tsat,sf ~ 4.9 is close to the saturation index value (about 4.3) of the soft 
component of RXTE X-ray spectrum from GRS 1915-1-105 (TS09). In ^ we argue that 
index saturations of the soft component at high m is expected for each given value of Ci/Cq. 
When Vma = 0.01(1 -|- Ci/Cq)c is low (or Ci/Cq ^5), the TL layer optical depth tq is high 
(see Eq. [6]) and consequently the plasma temperature of accretion fiow Tg is low (see Eq. 
[1]). In other words matter accumulation efficiency in the transition layer (or tq) is inversely 
proportional to magneto-acoustic velocity Vma there (see Eq. [2]). 
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ST09 and TS09 present strong observational and theoretical arguments that the satura- 
tion of index occurs at high m too. Specifically, ST09 argues that, in the converging flow, 
the spectral index, as an inverse of Comptonization Y— parameter, a ~ 1/Y = l/{7]Nsc), 
should saturate at high m because a linear increase of Ngc with m is compensated for by a 
decrease of 77 as 1/m. Here rj is average fractional energy per scattering. Our simulations 
show that the index saturation values of Tsat,h weakly depend on Ci/Cq ratio i.e. Tsat,h is in 
the range of 2.6-2.8 for all Ci/Co>5 (see Fig. E)- 

Particular for Ci/Cq = 5, the hard component of the simulated spectrum (above 80 
keV) is absent when rh increases from 0.1 to 0.5. (see Fig. [5]). However one can see that the 
spectral hard component arises when m increases above 1 (see also Fig. [6]). 

Moreover, for Ci/Cq = 5.5, and m^O.5 the flow electron temperature Tg decreases to 
about 3 keV (see Fig. [2]). In this case, we obtain that Th and Tgf saturate to ~ 2.7 and 
~ 4.3, respectively (see Fig. [9]). This situation is probably realized in the innermost part 
of GRS 1915+105 where there are similar values of ~ 2.7 and Ygf ~ 4.3 inferred from 
RXTE X-Taj spectrum (see TS09). 

One can conclude that for TL velocities Vma of order 0.05c and less, when the accu- 
mulation effect in the transition layer is very strong, the steep extended power law, as a 
spectral signature of the converging, is screened by the optically thick material of the tran- 
sition layer. In this case, the end of the spectral evolution from the low-hard state is the 
very soft state characterized by a strong blackbody like component (see the red and green 
simulated spectra in the upper panel of Fig. [5]). The soft spectral component is dominant 
at m > 0.5 and an extended hard tail is only seen in the spectrum for high values of m (see 
green spectral histogram in Fig. [5]). TS09 show that this type of spectral evolution is seen 
in RXTE observations of GRS 1915+105. 

However, for VAfA- 0.09c, the accumulation effect is not as strong in the TL and as a 
consequence the steep power law extended up to 300 keV is seen at higher level of count rate 
then that for Vma ~ 0.05c (compare the upper and lower panels in Fig. [5]). In this case, we 
see a classical spectral evolution from LHS to VHS through the intermediate state that is 



observed in many BH binaries [see e.g lRemillard fc McClintock I (j2006[ ). ST09]. 



We reproduce, using our Monte Carlo simulations, the observed spectral evolution in 
BH binaries when we assume the velocity proflle, approximately constant in the transition 
layer and followed by f ree-fall below 3-Rs - In fa ct, because the transition layer shrinks with 



mass accretion rate [see lTitarchuk et al. I (119981 )] one can assume the constancy of Vma there. 



Thus we proceed with our Monte-Carlo simulations of photon propagation in this kind of 
atmosphere. 
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However in the framework of our spectral evolution model, we do not obtain an extended 
power law with the index saturation of 2.1 that has been detected in HSS of Cyg X-1, GX 339- 
4, 4U 1543-47 (see ST09). This kind of the hard tail is reproduced in our simulations when we 
assume that there is a sub-Keplerian, relatively hot accretion flow with the same temperature 
as the transition layer, in addition to the disk accretion flow. In these simulations, the ratio 
of the sub-Keplerian to the Keplerian accretion rate is rhsub/mK = 50 [see Fig. [TOl) . This ide a 
of advection dominated accretion flow (ADAF), first introduced by iNarayan fc Yi I ( 1l994j ). 
is very popular in the astrophysical community. 



To illustr ate this, we use Cyg X-1, where it has been shown [see e.g. iPeterson I ( 119781 ). 



Kaper I ( 119981 )] that the X-ray source is powered mainly by accretion from the strong stellar 
wind of the supergiant star. There is a strong observational indication of the presence of 
sub-Keplerian flow in this source. The temperature of this flow T^^sub can be dictated by 
the photons Comptonized in the transition layer. Presumably when the disk mass accretion 
mK increases the sub-Keplerian mass accretion rate "frisub increases too. In order to have 
the relatively hot ADAF, one varies mx in the low value range. To show this, the variation 
of rriK from 0.05 to . 1 lea ds to a decrease of ~ T^^sub from 60 to 20 keV (see Fig. 
Laurent fc TitarchukI (119991 ) obtained that the photon index of the converging flow power-law 



tail Th saturates to ~ 2.1 when Tg, 



30 keV. 



Hence we obtain the observable saturation index of 2.1 assuming the existence of two 
component flows in the aforementioned BHs. 



6. Evolution of spectral efold energy Efoid with mass accretion rate m 

Dynamical (bulk) Comptonization spectra can be distinguished from thermal Comp- 
tonization spectra when there is information about the evolution of Efoid with rh. This 
information is inferred from X-ray observations of spectral transitions in BHCs, where the 
evolution patterns are drastically different between the two Comptonizations. In Figure [TTl 
we plot how the high energy efold Efoi^ energy varies with the mass accretion rate rh. In the 
bulk motion Comptonization case the efold energy Efoid decreases and then increases when 
mass accretion rate m^l (see Figure [TT]) . 

Note that in the thermal Comptonization spectra Efoid is proportional to the plasma 
temperature Tg {Efoid ~ 2kTg). Given that kT^ saturates to constant value of 2 keV for 
m^l (see Fig. El), -EfoM— ene rgy should also saturate to low values of Efoid ~ 3 — 4 keV [see 



Farinelli &: Titarchuk I ( l2010l ) for a similar situation in the neutron star case] 
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7. Conclusions 

We have studied, using Monte Carlo simulations, the evolution of X-ray energy spectra 
from the low-hard state to the high-soft state, and found that a hard spectrum should be 
observed at relatively low values of mass accretion rate independently of the velocity profile 
(see Fig. Further, for low velocities of the accretion flow in the transition layer (TL) 
V ~ 0.05c, there is plasma accumulation at high mass accretion rates that produces high 
optical depth and thus one can barely see the high energy tail. An emergent spectrum is 
characterized by a blackbody like shape of temperature on order of 1 keV and the relatively 
weak hard tail (see the green spectral histogram in the upper panel in Fig. [S]). This kind 
of spectra is observed in the very soft state of BH binaries. We also found that when the 
TL velocities are of order of (0.07-0. l)c at high mass accretion rate the extended power-law 
component is more pronounced in the emergent spectra (see the green spectral histogram in 
the lower panel of Fig. [5]). For all velocity values in TL, our results predict index saturation 
vs mass accretion rate (see Figs. UM^- Finally, we showed that the efold energy Efoid of the 
emergent spectra decreases and then increases as the mass accretion rate increases (see Fig. 
[TTl) . This is a specific property of BH spectra . This type of behavior of Efnid y s m (or vs 



photon index) has been recently discovered by iTitarchuk fc Shaposhnikov I ( l2010l ) in RXTE 
observations of XTE J1550-564. 

Note the decrease of the spectral efold energy vs mass accretion rate occurs when the 
high energy photons formed in the converging flow are not seen in the resulting spectrum. 
In this case, the spectral shape is only determined by thermal Comptonization taking place 
in the outer parts of Compton cloud. 

We appreciate the referee for careful reading and his/her valuable comments on the 
content of the manuscript. We also acknowledge the manuscript editing by Charlie Bradshaw. 
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Converging inflow Bulk comptonization 



Fig. 1. — A schematic view of the proposed geometry for thermal and bulk Comptonization 
regions in a source hosting a BH with power law-like emission at high energies. The bulk 
Comptonization plus thermal spectrum (bulk BMCl plus thermal BMC2) arises in the in- 
nermost part of the transition layer (TL), where the disk blackbody-like seed photons are 
(thermally and dynamically) Comptonized by the in-falling material. 
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Fig. 2. — Plasma temperature (solid lines) and optical depth (dash lines) of the transition 
layer vs m. 



1 



10 100 
Energy (KeV) 



1000 



Fig. 3. — The appearance of Bulk Motion Comptonization in the Monte-Carlo simulated 
spectrum. The spectrum is shown in units of keV (keV cm~^s~^keV~"^) with absolute values 
that are multiplied by an arbitrary constant. The dotted line corresponds to a simulation 
with an assumption that there is no bulk effect in the converging flow zone. The solid line 
results from the same simulation adding bulk velocity. The effect of bulk motion Comp- 
tonization is clearly seen above 20-30 keV. The parameters used are: mass accretion rate in 
Eddington units m = 5, flow electron temperature kT^ = 4.3 keV, Ci/Cq = 7 and q = 0.01. 
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Fig. 4. — An example of simulated bulk+tliermal Comptonization spectrum (red histogram) 
and the XSPEC model bmcl*highcutl +bmc2* highcut2. The spectrum is shown in units of 
keV (keV cm~^s~^keV~^) with absolute values that are multiplied by an arbitrary constant. 
The hard and soft components are shown by dash green and dash red curves, respectively. 
The resulting spectrum is shown by a solid green curve. In the simulations, the mass accretion 
rate in Eddington units is m = 5, Ci/Cq = 5, g = 0.01 and the flow electron temperature is 
kTe = 2.4 keV. 
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Fig. 5. — Spectral evolution as a function of mass accretion rate: upper panel: Ci/Cq = 5, 
bottom panel: for Ci/Cq = 10. Black, blue, red and green histograms correspond to m = 
0.05, 0.1, 0.5 and 5, respectively. 
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Fig. 6. — Index vs m for two BMC components: red points are for the soft component and 
blue points are for the hard BMC component Cl/C0=5. 
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Fig. 7. — Same as that in Figure [6] but for Cl/C0=8. 
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Fig. 8. — Same as that in Figure [6] but for Cl/C0=10. 
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Fig. 9. — Saturation values of index as a function of Cl/CO. 
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Fig. 10. — Index vs m in the presence of sub-Keplerian flow above the disk. The accretion 
flow has the same electron temperature as the transition layer one. The ratio of rhsub/^K = 
50. Total rhtot = rfiK + rhsub 
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Fig. 11. — Efoid energy vs m: Ci/Cq = 5, 8, 10 upper, middle and bottom panels respec- 
tively. 



